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We apply a wilsonian renormalization group approach to the system of electrons in a two- 
(f) ' dimensional square lattice interacting near the saddle-points of the band, when the correlations at 

momentum Q = (tv, n) prevail in the system. The detailed consideration of the spin degrees of 
freedom allows to discern the way in which the SU(2) spin invariance is preserved in the renormal- 
ization process. Regarding the spin correlations, we find two different universality classes which 
correspond, in the context of the extended Hubbard model, to having the bare on-site interaction U 
repulsive or attractive. The first class is characterized by a spin instability which develops through 
the condensation of particle-hole pairs with momentum Q, with the disappearance of the Fermi line 
in the neighborhood of the saddle-points. Within that class, the attractive or repulsive character 
of the nearest-neighbor interaction V dictates whether there is or not a d-wave superconducting 
instability in the system. For the Hubbard model with just on-site interaction, we show that some 
of the irrelevant operators are able to trigger the superconducting instability. The naturalness of 
the competing instabilities is guaranteed by the existence of a range of doping levels in which the 
chemical potential of the open system is renormalized to the level of the saddle-points. We incor- 
porate this effect to obtain the phase diagram as a function of the bare chemical potential, which 
displays a point of optimal doping separating the regions of superconductivity and spin instability. 
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I. INTRODUCTION 



During the last years there has been much effort devoted to the study of strongly correlated electron systems. The 
I ' interest has been maintained by the behavior displayed by the high-T c copper-oxide compounds since the discovery of 
, their superconductivity 15 years ago [1]. There are a number of features exhibited by these materials that do not fit 
into the conventional theoretical frameworks. The normal state of the cuprates shows for instance unusual transport 
properties and, more strikingly, a pseudogap phase in which part of the density of states is lost at the Fermi level 
while the system remains conducting. It seems that a new paradigm is needed to describe these materials, in the 
same way as the Fermi liquid picture accounts for the behavior of conventional metals. 

From the theoretical point of view, progress has been made during the past decade in understanding the foundations 
of Landau's Fermi liquid theory and, consequently, the possible deviations that may open the way to a new kind of 
metallic behavior [2-4]. The most powerful method used in this task has been the renormalization group (RG) 
approach developed for interacting fermion systems [2]. We have learned from it that the Fermi liquid picture is 
a very robust description of the metallic state. There are only a few perturbations that may destabilize the Fermi 
I ■ liquid, favoring the formation of states with different types of symmetry breaking. The Fermi liquid represents itself 
a universality class in which any electron system falls at dimension D > 2, unless the interaction is sufficiently 
long-ranged [5—11] or the Fermi surface develops singular points [12]. 

Soon after the discovery of the high-T c superconductivity, it was proposed that the presence of nonlinear dispersion 
near the Fermi line of the copper-oxide layers could be at the origin of the unconventional behavior [13,14]. The 
fermion systems in a two-dimensional (2D) square lattice have necessarily saddle-points in their band dispersion, 
which give rise to Van Hove singularities where the density of states diverges logarithmically. In the most common 
instances, the two inequivalent saddle-points lie at the boundary of the Brillouin Zone, and their hybridization has 
been proposed to explain the existence of a ci-wave order parameter in the superconducting phase [15-19], as observed 
experimentally. Further investigations have shown that the unconventional transport properties in the normal state 
may be accounted for by the proximity of the Fermi level to the Van Hove singularity (VHS) in the copper-oxide 
layers [20-24]. 

A careful examination of the kinematics near the saddle-points has shown indeed that a superconducting instability 
with d-wave order parameter arises in the t—t' Hubbard model with bare repulsive interaction [25,26]. The mechanism 
at work is of the same kind described by Kohn and Luttinger as giving rise to a p-wave pairing instability in the 
three-dimensional Fermi liquid [27,28], but adapted now to the 2D model with saddle-points near the Fermi line. 
Other studies have considered in detail the influence of the entire Fermi line in the development of the instabilities 
of the system [29,30]. They have given further support to the picture of a competition between a spin-density- wave 
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instability and a pairing instability with d-wave order parameter in the t — t' Hubbard model with the Fermi level 
at the VHS. More recently, a refined renormalization program has been implemented in Ref. [31] by trying to handle 
the momentum dependence of the vertex functions in the scaling procedure, what has confirmed the appearance of 
different phases with symmetry breaking in the spin and the charge sector. 

Despite all the results obtained in the system of electrons near the VHS, there are still important obstacles precluding 
a precise description of the effective theory at low energies. From a technical point of view, the source of the problem 
is the appearance of infrared singularities in the RG approach after accomplishing the renormalization of the leading 
logarithm. Some vertex functions, like the four-point interaction with vanishing total incoming momentum at the 
one-loop level or the electron self-energy at the two-loop level, get log 2 (A) corrections in terms of the energy cutoff 
A. After applying the standard RG program, the renormalized quantities still contain factors of the form log(A). 
This fact questions the predictability of the theory since the argument of the logarithm has a hidden energy scale, 
which sets the strength of the corrections. From a formal point of view, the theory becomes nonrenormalizable in 
the standard RG approach, since the energy cutoff is not the only dimensionful variable that appears in the scaling 
process. 

The problem of the rcnormalizability of the theory can be best handled by adopting a wilsonian RG approach, 
in which only the high-energy modes that live at the cutoff A are integrated out at each RG step. In the present 
paper we follow Shankar's RG program for interacting fermion systems [2], which has the advantage of decoupling 
the renormalization of the BCS channel (with vanishing momentum of the colliding particles) from that of the rest of 
the channels at the one-loop level. 

Moreover, the important feature of the wilsonian approach is that it allows to set free the chemical potential, so 
that it can readjust itself at each step of integration. The issue of the renormalization of the chemical potential has 
been discussed in Ref. [2] in the context of Fermi liquid theory, and it reaches great significance when considering the 
system of electrons near the VHS. The chemical potential cannot be fixed at the singularity from the start, since it is 
actually the scale needed to regularize the infrared singularities that appear in the standard RG procedure. On the 
other hand, the final location of the chemical potential relative to the VHS is not arbitrary, since it is a dynamical 
quantity that scales in a predictable way upon renormalization. 

We remark that the renormalization devised in the paper assumes a constant value of the bare chemical potential, 
instead of a constant particle number of the system. That is, we describe a situation appropriate for an electron 
system in contact with a charge reservoir, which sets the nominal value of the ensemble. The renormalization 
accounts for the reduction suffered by the effective chemical potential inside the electron system due to the repulsive 
interaction. This description of the electron system at constant nominal chemical potential is most appropriate when 
dealing with the Cu-0 layers of the cuprate superconductors, since it provides a realization of the contact of the 2D 
layers with the charge reservoir. The conclusion is that a variation in the external chemical potential does not have 
always a linear correspondence with the variation of the final renormalized value of /i, which is identified with the 
Fermi energy of the electron system. 

The renormalization of the chemical potential makes possible to address the question of the naturalness of the 
picture in which the Fermi level is fine-tuned to the VHS. The strength of the predicted instabilities depends crucially 
on the proximity of the Fermi energy to the singularity. This has been the main criticism to the proposals claiming 
that the features of the copper-oxide materials could be related to the properties of electrons interacting near a VHS. 
We will show that the chemical potential is renormalized towards the VHS in a certain range of filling levels, in such 
a way that it may become pinned to the singularity in the low-energy theory. This fact was already anticipated in 
Refs. [32], [33] and [16], and it has been used to cure the infrared singularities of the electron self-energy in Ref. [21]. 
In the present paper, we will take into account such an effect to determine in a predictable way the strength of the 
pairing instability in the system, as a function of the different values of the bare chemical potential. 

In the next section we describe the system to which our analysis applies. In Section III we classify the different 
renormalized vertices that arise by explicit consideration of the spin degrees of freedom. The universality classes of 
the system are obtained in Section IV, where we also show the way in which the SU(2) spin invariance is preserved 
along the RG flow. Section V is devoted to establish the properties of the spin instability of the system, while Section 
VI analyzes the renormalization of the chemical potential to determine the region of the phase diagram in which the 
superconducting instability prevails. Finally, the last section is devoted to draw the main conclusions of this work. 

II. THE MODEL 

We take as starting point of our analysis a system of interacting electrons in the 2D square lattice with nearest- 
neighbor hopping t and next-to-nearest-neighbor hopping t' . The band dispersion of the model is given by 

e(k) = — 2t(cos(k x ) + cos(k y j) + At'cos(k x )cos(k y ) (1) 
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where we have set the lattice spacing equal to one. Some of the energy contour lines are shown in Fig. 1. The 
dispersion has two inequivalent saddle-points A and B at the boundary of the Brillouin Zone. In their neighborhood, 
the energy of the one-particle states can be approximated by the quadratic form 

£A,u(k) ^T(tT2t')k 2 x ±(t±2t')k 2 y (2) 

where the momenta k x and k v measure now small deviations from A and B. 




As a consequence of the nonlinear character of the dispersion, the density of states n(e) diverges logarithmically at 
the level of the saddle-points 

n(s)^clog(t/\e\)/(4n 2 t) (3) 



with c = 1/^/1 — 4(i'/i) 2 . This implies that, when the Fermi level is close to the VHS, most part of the low-energy 
states are concentrated in the neighborhood of the two saddle-points. In order to apply the RG approach, we may 
take two patches where the quadratic approximation (2) holds around the saddle-points. Higher-order corrections 
to the expression (2) are irrelevant under the scaling that makes the action of the model a fixed-point of the RG 
transformations, as we see in what follows. 

We consider then a model whose action at the classical level is 

S = ^ f dtd 2 p (z*+ ff (p)d t * aff (p) - ( £o (p) - mo) *+ (P)*aa(p)) 

+ y > * d Vrf 2 P2rf 2 P3rfV^(Pl,P2,P3,P4)*^(Pl)*^(P2)*c CT KP4)*d ( r(P3)'5(Pl +P2 -P3 -P4) (4) 

a.b.c.d 

where the indices a, 6, c, d run over the two patches around A and B. 

The scaling transformation that leaves invariant the kinetic term of the action is 

d t - sd t (5) 
p-aVSp (6) 

*a.(p) - .S- 1/2 * aCT ( P ) (7) 

It is easily checked that, with the transformation (5)-(7), the interaction term in the action (4) is also scale invariant 
for a constant value of the potential U(pi, P2, P3, P4). If this is not constant, provided that it is a regular function 
of the arguments we can resort to an expansion in powers of the momenta. Only the constant term is significant, 
since the rest of higher-order terms fade away upon scaling to the low-energy limit s — > 0. This means that we meet 
the first requirement to apply the RG program, that is to have a model which converges to a fixed-point under RG 
transformations at the classical level. 

In the above scaling, we already find the first deviation in the RG program with respect to the analysis of Fermi 
liquid theory. In the case of a model with circular Fermi line, the interaction term is scale invariant only for very 
special kinematics of the scattering processes [2]. In our model, we have seen that no constraint is needed on the 
four momenta involved in the interaction at the classical level. It is only after taking into account virtual processes 
that the interactions will start to grow large under scaling for some particular choices of the kinematics. This will 
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single out a number of so-called marginally relevant channels among all the scattering processes, recovering then the 
similitude with the analysis of Fermi liquid theory at the quantum level. 

The two-patch RG analysis of the t — t' Hubbard model has proven to give the dominant instabilities of the system 
with the Fermi level at the VHS. For t' > 0.276 t, a ferromagnetic phase has been found below a certain critical 
frequency [25,26,34,35], in agreement with the results obtained from Monte Carlo calculations [36]. In this paper we 
will be interested in the regime with t' < 0.276 t, where the competition between a spin instability and a pairing 
instability arises, making the model more appropriate for the comparison with the phenomenology of the cuprates. 



III. WILSONIAN RENORMALIZATION GROUP 



In what follows we apply a wilsonian RG approach to obtain the low-energy effective theory of the system. We 
proceed by progressive integration of the modes in two thin shells of width dA at distance A in energy below and above 
the Fermi level, as depicted in Fig. 2. For the time being, we will assume that the Fermi level is located precisely 
at the VHS, unless otherwise stated. This is crucial to obtain a significant renormalization in any of the interaction 
channels, and later on we will comment on the naturalness of this situation. 




FIG. 2. Picture of the density of states n(e) and of the renormalization of the chemical potential fj, by integration of states 
at the energy cutoff A. 



The vertex functions may become relevant, that is increasingly large at low energies, only for very definite choices 
of the kinematics. Focusing on the four-point interaction vertex, this is renormalized by a quantity of order dA at 
each RG step only when the momentum transfer along a pair of external lines is either or Q = (tt,tt), or when 
the total momentum of the incoming modes vanishes (BCS channel). In the present work we deal with the latter 
two instances, since the first corresponds to the case of forward-scattering interactions, which are subdominant in the 
range t' < 0.276 t that we are considering. In this regime, the divergences at vanishing momentum-transfer are related 
to charge instabilities of the system, which have been treated in detail elsewhere [37]. We will see that divergences 
in the channel with momentum transfer Q give rise to a spin instability, which competes with the superconducting 
instability in the BCS channel in the model with a bare on-site repulsive interaction. 

The different kinematics which may appear in the BCS channel are listed in Fig. 3. We allow for the possibility of 
Umklapp processes in which the incoming modes scatter from one of the saddle points to the other. 





FIG 



(a) (b) 

3. BCS vertices that undergo renormalization by particle-particle diagrams. 



The solid and dashed lines stand for modes 



in the neighborhood of the two different saddle points. 
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The different kinematical possibilities that arise in the channel with momentum transfer Q are classified in Figs. 
4 and 5. The first includes the interactions in which the incoming modes are at different saddle points, while the 
latter contains the Umklapp processes. The other important distinction is between direct (D) and exchange (E) 
interactions. Direct processes are those in which the momentum transfer Q is taken by the same scattered fermion 
line, while in a exchange process the momentum transfer takes place between two different fermion lines connected 
only by the interaction. 



k+Q t 



P-Q t 




k+Q t 






FIG. 4. Direct and exchange vertices that undergo renormalization by particle-hole diagrams. 



k+Q t 




k+Q t 



p+Q j 




p+Q t 




p+Q t 




FIG. 5. Umklapp vertices that undergo renormalization by particle-hole diagrams. 



The interaction vertices depicted in Figs. 3-5 are all renormalized upon reduction of the cutoff A. This can be 
traced back to the divergent behavior of the different susceptibilities of the model. By integration of the high-energy 
modes in the shells of width dA, the particle-hole susceptibility at momentum Q gets a contribution 



(8) 



where d = log ^1 + ^/1— 4(t'/t) 2 ^j /(2t'/t) [38]. In the same fashion, the contribution to the particle-particle 
susceptibility at zero total momentum is 
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4irH 



log(A)dA/A 



(9) 



In the latter case, the result of the differential integration diverges logarithmically in the limit A — > 0. This has 
been a source of problems in the usual RG analyses of the model. The definition of the argument in the logarithm 
needs an additional scale, while a proper RG scaling requires that the energy is the only dimcnsionful variable in the 
problem. It has to be realized that the coefficient at the right-hand-side of Eq. (9) represents actually the density of 
states. This has to be born in mind for the correct implementation of the RG approach, as we will discuss later. 

Let us deal first with the renormalization of the vertices with BCS kinematics in Fig. 3. At the one-loop level, the 
vertices Vj and Vjj get corrections of order dA/A from the diagrams shown in Fig. 6. It is important to realize that 
these are the only diagrams to be taken into account to first order in dA. There are also corrections from particle-hole 
diagrams but, as long as the momentum that goes into the particle-hole loop is not precisely zero or Q, these terms 
are of order (dA) 2 and therefore irrelevant in the low-energy limit, as shown graphically in Fig. 7. 





P P P P 

(a) (b) 





P P P P 

(c) (d) 

FIG. 6. Particle-particle diagrams renormalizing the BCS vertices at the one-loop level. 









/ / 


A — ^ 















FIG. 7. Picture of the high-energy shells of width dA at a given saddle-point. The dark regions represent the contribution 
to a particle-hole diagram when q is the total incoming momentum. 



The BCS vertices mix between themselves alone at the one-loop level, and the situation is similar in that respect 
to the general analysis of the 2D Fermi liquid [2]. The degree of renormalization depends on the density of states 
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n(e) at the shells integrated out. For later use, we consider at this point the most general case in which the chemical 
potential \i does not coincide from the start with the level of the VHS. The differential RG equations take then the 
form 

A^=cn( M -A) (V? + V*) (10) 
A-^- = 2c n(p - A) VjVu (11) 

These equations were considered in Ref. [25], and they also appear as the leading order in the RG approach of Rcf. 
[31]- 

We consider next the renormalization of the vertices -Eqj_ arid Eu±, which have also the property that they mix 
only between themselves in the one-loop corrections linear in dA. These have been represented in Fig. 8. It can be 
checked that any other diagrams give irrelevant contributions of order (<iA) 2 , because they involve either a particle- 
hole susceptibility at momentum different from Q or a particle-particle susceptibility with total momentum different 
from zero. In the latter case, for instance, it is shown in Fig. 9 that the number of intermediate states produced by 
integration of high-energy modes is quadratic, instead of linear in dA. 

p-Q t k+Ql p-Q t k+QI 

\ /' \ / 



kt Pi kt Pi 

(a) (b) 

p+Qt k+Ql p+Qt k+Ql 

\ / \ / 




kt Pi kt Pi 

(c) (d) 

FIG. 8. Particle-hole diagrams renormalizing the vertices Eq± and Eu± at the one-loop level. 




FIG. 9. Same scheme as in Fig. 7. The dark regions represent the contribution to a particle-particle diagram when q is the 
total incoming momentum. 
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The differential RG equations for the pair of vertices read 

A^ = -c'(E 2 Q± +E* u± )/(4nH) (12) 
A ^lr = - c ' E ^ E uJ(^ 2 t) (13) 

These equations were obtained in Ref. [25], where the names Ui n t er and V um k were used instead of Eq± and Ejj±_ 
introduced in the present paper. The same equations also arise at the dominant level in the functional renormalization 
of Ref. [31]. 

We now turn to the rest of the vertices, Dqm , Dq±, Eq\\ , Djj\\ , Djj± and Eu\\ , which renormalize among themselves 
at the one-loop level. It is clear that the vertices Dqy and Eq\\ cannot be distinguished from each other just by looking 
at the external legs. The same applies to Du\\ and Ejj\\. At the one-loop level, one can still discern whether the 
momentum transfer Q takes place along the same scattered fermion line or not. However, the different corrections 
have to organize so that the above pairs of vertices enter in the combinations Dqu — Dqu and Dj/|| — E^, which are 
the quantities that make physical sense. In that respect, the situation is similar to what happens with the couplings 
<7iH and §2\\ in the one-dimensional electron systems [39]. 

The one-loop renormalization of the vertices provides an explicit proof of the above statement. The vertex Dqu 
gets linear corrections in dA from the diagrams shown in Fig. 10, while Eq\\ is rcnormalized by the diagrams shown 
in Fig. 11. Their RG equations read then 



3D, 

A 



' { D Q\\ + D Q^ + D u\\ + D u± - 2D Q\\ E Q\\ - 2^11^1,) /(^ 2 t) (14) 



OA 

A ^ = - c '{ E *\\ +E h)/^ 2t ) ( 15 ) 

(16) 



These two equations can be combined to be written in terms of the physical vertex, 

8 (D oll - Er 



a - v^QII ~ ^Qll 
dA 



(Dqh-Dqh) +(D mi -E u{l ) +D 2 Q± + D 2 U± /(AttH) (17) 



k+Qt p-Qt k+Qt p-Qt 




k+Qt p-Ql k+Qt p-Qt 




(c) (d) 



FIG. 10. Particle-hole diagrams renormalizing the vertex Dqu at the one-loop level. 
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P-Qt 



k+Qt p-Q t 



k+Qt 





FIG. 11. Particle-hole diagrams renormalizing the vertex Eq« at the one-loop level. 

The RG equations for the remaining vertices also depend on the combinations Dq|| — Eq\\ and Du\\ — E v \\. In the 
case of Dq±, we have 



A 



= c ' [v°Q|| - ^Qll) D ^ + ( D m\ E U{1 ) D U± ] /(2nH) 



Finally, the RG equations for £>c/||, Djj± and Ejj\\ take the form 

A — ^ = c' (£>q|| - D m E m - A/||£ Q || + D Q± D U± ) /(2n 2 t) 



A- 



<9d, 



u\\ 



dA 



-c'E Qll E ull /(2n 2 t) 



A 



As a final check, the equation for Du\\ — E v \\ turns out to depend on the physical combination of couplings 
d{D m -E v {) 



A- 



<9A 



= c' [(D Ql{ - E Q ||) (D^h - E m ) + D Q ±D UA ] /(2ttH) 



(18) 

(19) 
(20) 
(21) 

(22) 



IV. UNIVERSALITY CLASSES 



We discuss now the universality classes in which the system may fall regarding the spin correlations. We will 
focus on the analysis of bare repulsive interactions, that is where the competition between spin and superconducting 
instabilities arises. We will see that our RG scheme is able to preserve the spin-rotational invariance of models whose 
bare interactions have such a symmetry. This provides another nontrivial check of our RG approach, as our framework 
offers the possibility to analyze the scaling of interactions with and without the SU(2) spin symmetry. 

The interactions of physical interest have the property that Dq|| — Dq|| = D{/|| — Ejj\\ and Dq± = Djjj_- These 
conditions are maintained along the RG flow if they are satisfied by the bare couplings. Thus, it is useful to work 
with the set of couplings 



Df =Dr 



E m ±D L 



E, 



Df = D Q± ± D u± 



From the results of the preceding section, these new couplings satisfy the equations 



A 



A 



8A 
dA 



(Df)* + (Dff\/(4nH) 
= c'D±D±/(2ir 2 t) 



(23) 
(24) 



(25) 
(26) 



The universality classes of the system can be obtained from the integrals of Eqs. (25) and (26). We stick to the case 
in which D^ = DJ = 0. The flow for the couplings Djj~ and D~\_ is represented in Fig. 12. Focusing on interactions 

that are repulsive at the initial stage of the RG, that is Djf > and D~\_ > 0, we observe two possible behaviors of 
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the renormalized couplings. In the case in which the bare couplings satisfy > D^, the flow is bounded and it 
converges monotonically to the origin of the space of couplings. If we start otherwise from a point with < Z)+, 
the flow becomes unstable and it approaches a regime in which Dt — > — oo and D~\_ — » +00. 




FIG. 12. Flow of the renormalized couplings in the (Dt ,D~J_) plane 



The regions with stable and unstable flow correspond to respective universality classes, which imply quite different 
physical properties. Let us focus, for instance, on the extended Hubbard model with on-site interaction U and 
interaction V between nearest-neighbor sites. The appropriate bare values for the couplings in Figs. 4 and 5 are 







= U 


-w 


Eq\\ 


= Eq± 


= U 


+ aV 




= D U± 


= U 


-W 


Em\ 


— Eu± 


= u 


-0V 



(27) 
(28) 
(29) 
(30) 



with < a, [3 < 4. We have for the initial values of the flow Djj" = —(8 + a — (3)V and D+ = 2U - 8V. With the 
physically sensible choice a = (3, we see that the attractive or repulsive character of the on-site interaction dictates 
whether the RG flow is bounded or not in the upper half-plane of Fig. 12. 

The fact that the flow is not bounded for U > points to the development of some instability in the system. The 
divergence of the renormalized couplings represents the failure to describe the model in terms of the original fermion 
variables. The underlying physical effect is the condensation of boson degrees of freedom, as we will show in the next 
section. The preservation of the spin-rotational invariance at each step of the RG process helps to clarify the physical 
interpretation of the instability and to discern the issue of the spontaneous breakdown of the symmetry. 

We pay attention then to the way in which the SU(2) spin symmetry is preserved in our RG framework. This 
can be analyzed by looking at the response functions for the different components of the spin operator. Since the 
renormalized interactions grow large at momentum transfer Q = (it, it), we focus on the correlations of the operator 



Sj(Q) 



^^(k + Q^'vMk) 



J =x,y,z 



(31) 



The scaling properties of the response functions can be studied in the same fashion as for the interacting one- 
dimensional fermion systems [40]. The response function R z (u>) for the 5 Z (Q) operator, for instance, is renormalized 
by the diagrams shown in Fig. 13. After taking the derivative with respect to the cutoff and imposing the self- 
consistency of the diagrammatic expansion, we obtain 



dR z 
OA 



2c' 1 c' 



Eq\\ + Djj\\ - Eu\\ - Dq± - Du_i) —R 



(32) 
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P+Qt 



(a) 



q+Qt 




P+Qt 



(b) 



q+Qt 




P+Qt 



(c) 



q+Q| 




P+Qt 



(d) 



q+Q| 



P+Q(X 



qt 



q+Qt 



P+Qt 




(e) <0 

FIG. 13. First-order contributions to the correlator of the S z operator. 

The response functions R x {lo) and R y (co) for the other two components of the spin operator are both renormalized 
by the diagrams shown in Fig. 14. Following the same procedure as for R z (ui), we obtain 



dRx 2C ' 1 C ' (F +F 



and a completely similar equation for R v (u>). 



(33) 



P+QiX 



q+Qj 




(a) (b) 
FIG. 14. First-order contributions to the correlators of the S x and S y operators. 

The response functions R x (ui), R y (to) and R z (ui) can be made exactly equal if the equation 



Dr 



Eq\\ + D v \\ 



E, 



u\\ 



D, 



Q-L 



Du± — —E{ 



Q_L 



(34) 



is satisfied all along the flow. From Eqs. (12), (13), (25), and (26), we observe that this is automatically fulfilled 
when the condition is imposed for the initial values of the couplings. In the case of the extended Hubbard model, we 
have indeed for the bare couplings in Eqs. (27)-(30) 



Dq± + D u± - D Q \\ +E Q \\- D m +E V \\= E Q± + E u± = 2U + (a - (3)V 



(35) 



The condition is actually satisfied by the couplings of any hamiltonian that is invariant under rotations. We show in 
this way that the SU(2) spin symmetry can be preserved at each point of the RG flow of the couplings, so that the 
low-energy effective action keeps the invariance of the bare hamiltonian. 



V. SPIN INSTABILITY 



We proceed to determine the physical properties of the universality class corresponding to the unstable flow in 
the upper half-plane of Fig. 12. The divergence of the renormalized couplings — Dt and Eq± + Ejj± results 
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in the divergence of the response functions R x , R y and R z at a certain value of their argument. This points at the 
development of an instability in the spin sector at the corresponding value of the energy measured from the Fermi 
level. 

The divergence of the response functions implies the existence of a pole at a given frequency u) c . From the solution 
to Eqs. (12), (13), (25), and (26), the value of the pole is given by 

1 - (D+(A ) - D+(Ao))x P h(Q^ c ) = (36) 

where D^(Ao) and Djj~(Ao) are the initial values of the couplings. As long as the susceptibility Xph at momentum 
Q diverges logarithmically in the low- frequency limit, it is clear that the above condition is satisfied no matter how 
small the initial value of the coupling — may be. 

It is important to bear in mind that the susceptibility Xph at momentum Q has a finite imaginary part, which is 
essential to discern the nature of the ground state of the system. The imaginary part is computed in the Appendix, 
and it turns out to be c'/(8irt). The equation (36) can be written then in the form 

1 - (D+(A ) - D+(A ))^ log(zAoK) = (37) 

which shows that the pole occurs for a pure imaginary value oj c = i\ui c \. 

The appearance of a pole in the correlator of a boson operator for a pure imaginary frequency corresponds to a 
phenomenon of condensation, in the same fashion as it happens in the case of a pairing instability [41]. In the present 
instance, the boson-like object is the spin operator at momentum Q defined in Eq. (31). The fact that the pole arises 
at a value i\w c \ means that the instability pertains actually to the theory posed at finite temperature, and that there 
is a transition to a condensed phase at a temperature of the order of magnitude given by \u c \. 

In our case, the boson operator that acquires a nonvanishing mean value due to the spin instability is the vector 
/ d 2 kdu)^+{k)(T aa 1v(k + Q). This has important consequences, since the diagrammatic approach has to be rebuilt 
below the point of the transition, in the same way as in the case of a pairing instability [42] . 

Let us focus on the Hubbard model, i. e. on a model with interaction between currents with opposite spin 
projections. To fix ideas, suppose that the vector S gets the nonzero mean value pointing in the x direction. Then, 
there arc two different kinds of one-particle propagators, since the presence of the condensate leads to the consideration 
of correlators of the type (^^(k, w^sj^k, u;)), as well as of the usual propagators for well-defined spin projection 
near each of the saddle-points. To include all the different possibilities, wc define the propagator G a(Tj b (7 '(k, w), with 
indices a, b labelling the saddle-points and a, a' labelling the spin projections: 

G aaM ,(k,u) = i<*+ (k.wjtf^k.w)) (38) 

The Schwinger-Dyson equations for the one-particle propagators take the form shown graphically in Fig. 15, where 
the insertion of the wavy line represents the factor 

C/y*d 2 ^(*+ T (k)* Bi (k + Q)) = A (39) 
We have, for instance, the closed set of equations 

G A ]M = G %,A\ + G A].At AG Bl,A] (40) 

Gb U m=G% bi A*G a ^ m (41) 

where the superindex denotes the corresponding propagator before the introduction of the condensate. Eqs. (40) 
and (41) can be combined to give an equation for G A -\,A-\, which reads 

G A ],A] = G A ],A1 + G A1,A1 AG B\,Bl A * G AhA1 (42) 
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FIG. 15. Self-consistent equations for the dressed propagators in the particle-hole condensate, in terms of the undressed 
propagators at the two inequivalent saddle-points. 

The solution to Eq. (42) takes the form 

Ga\.a\ (k, to) = G aIai(&u) (43) 

l-G^ ijlT (k,u;)|ApGg >Bi (k, W ) 



in terms of the propagators at the two different saddle-points 

g<°> T (k, W ) = ■ ■ - — ; - — - ( j -i i 



'AiM**")- w -e A (k)+ie sgnH 
,(o) , 1 



Go B . (k, w) = — (45) 

The important point is to determine the pole structure of the propagator (43). Its frequency dependence can be 
expressed in the form 

Gaim (k, to) = ^ - ^ (k) - j£ ggn(w)) ^ - ^ (k) - . £ ggn(w)) - |A|2 (46) 
«( k ) 2 , ^( k ) 2 (47) 



w — e M (k) + ie sgn(w) w — e„(k) + it sgn(w) 
with appropriate weights u(k) 2 ,w(k) 2 , and e„(k),e„(k) being the roots of the denominator in Eq. (46) 

s u ,v (k) - (s A (k) + eb (k) ± V(£A(k)-e B (k)) 2 +4|A| 2 ) /2 (48) 

From the physical point of view, the most important feature is the appearance of a gap in the quasiparticle spectrum 
near the saddle-points. This can be checked by determining the shape of the Fermi line, which is given by setting 
either e„(k) = or e„(k) = 0. Both conditions lead to the equation 

e^(k) £B (k) - |A| 2 = (49) 

By recalling that £A(k) = + t+ky and £s(k) = i+fc 2 — t-ky, we end up with the equation satisfied by the 

points of the Fermi line 

(t_fc 2 -t + fc 2 )(t + fc 2 -t_fc 2 ) + |A| 2 =0 (50) 
Solving Eq. (50) for the variable fc 2 , for instance, we find that there is a solution only for values of A; 2 such that 

(4-t 2 _) 2 fc^-4|A| 2 M- >0 (51) 

Reminding that t± t ± 2t', this condition implies that, for small values of t', there is a gap in the spectrum of 
quasiparticles in the range 
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M% Z |A| 



(52) 



We see therefore that the gap opens u p in the neighborhood of the saddle-points. The size of the part of the Fermi 
line destroyed is bounded by y/\A\/t', in units of the inverse lattice spacing. 

The formation of the quasiparticle gap has its origin in the hybridization of modes at different saddle-points, as a 
consequence of the enhanced scattering with momentum transfer exactly equal to Q. Quite remarkably, this is an 
effect that can be studied in the weak coupling regime of the model, and the gap appears for arbitrarily small strength 
U of the interaction. From the technical point of view, the discussion carried out in this section parallels the treatment 
of the one-particle Green functions in the usual description of the superconducting instability [42]. However, it is clear 
that the physical setting is quite different. In the present situation, the condensate is made of particle-hole pairs with 
a nonvanishing average projection of the spin. The fact that a macroscopic number of these pairs has been formed is 
what forces the quasiparticles to live out of the range already excited by the condensate. 

An important issue concerns the spontaneous breakdown of the spin-rotational symmetry in the condensate. Let us 
consider the model at zero temperature regarding this matter. It is clear that the nonvanishing average spin cannot 
have in principle any preferred direction in space. Recalling our definition in Eq. (39), a real value of A implies that 
the spin of the condensate points in the x direction, since 

J d 2 fc^(*+ T (k)*Bi(k + Q)> + J d 2 fcdw(*+ i (k + Q)*At(k)> + A «-» B = 2(A + A*)/U (53) 

A purely imaginary value of A implies otherwise that the spin of the condensate lies in the y direction. Finally, it 
may also be that the nonvanishing mean value is realized for the z component of the spin 

J d 2 fcrfo;(*+ T (k)* BT (k + Q)) - J d 2 fcdo;(*+ i (k)* i ji(k + Q))+A^ J B^0 (54) 

In the ground state of the model at zero temperature, the spin of the condensate has to point in a definite direction 
and the SU(2) rotational symmetry is spontaneously broken. As a consequence, two Goldstone bosons arise in the 
spectrum, which correspond to the spin waves that propagate on top of the particle-hole condensate. These are the 
gapless excitations of the model, together with the quasiparticle excitations that exist sufficiently far away from the 
saddle-points. 



VI. SUPERCONDUCTING INSTABILITY 



We now turn to the instability that arises from the divergent flow of Eqs. (10) and (11). The integral of these 
equations depends on the position of the chemical potential with respect to the VHS. For this reason, it is crucial to 
know how fi depends on the cutoff A as this is progressively lowered. 

The issue of the renormalization of the chemical potential has to be treated necessarily in the framework of the 
wilsonian RG approach. As the high-energy modes are integrated out at the scale A, shifts its position by a quantity 
propotional to dA. At the same time, it is the chemical potential which sets the level to measure the energy cutoff, 
as shown graphically in Fig. 2. The outcome is that fi adjusts itself at each step of the RG process, until the point 
in which the cutoff A is lowered down to the final chemical potential. 

At the computational level, the shift of \i is given by the frequency and momentum-independent part of the electron 
self-energy, with intermediate states taken from the high-energy modes being integrated. The renormalization is 
proportional to the charge of the occupied states in the lower slice of width dA, which couples through the forward- 
scattering vertex F in the usual Hartree and exchange diagrams. The RG equation for the chemical potential reads 

^=F(/i-A)n( M -A) (55) 

The perturbative approach is further improved by incorporating the renormalization of the F vertex, which bears a 
well-known dependence on the energy scale measured from the VHS [26,37] 

F(e) w F /(l - F log(|e|)/(4^)) (56) 

When the density of states n(e) is a smooth function of the energy, the integration of high-energy modes produces 
a steady downward flow of \x. The physical interpretation of this effect corresponds to the upward displacement of 
the one-particle levels due to the repulsive electronic interaction. In the neighborhood of the VHS, the dynamics of 
\i becomes highly nonlinear given the singular behavior of the density of states in Eq. (55). It turns out that, in 
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a certain range of initial values, the chemical potential is renormalized down to the VHS and precisely pinned to it 
in the low-energy regime. As stated in the Introduction, this result pertains to a statistical description in terms of 
the grand canonical ensemble. The physical picture is appropriate then for an open system in contact with a charge 
reservoir, which sets the bare value [io of the chemical potential. 

In order to evaluate the influence of the VHS on the renormalization of the chemical potential, we have solved Eq. 
(55) with the approximate density of states 



n{e) = c\og{t/\e\)/{ATT 2 t) for |e| < 0.5i (57) 
const. for \e\ > 0.5i (58) 



This expression has the correct normalization for the logarithmic singularity in the 2D square lattice. The behavior 
of the integrals of Eq. (55) with such a density of states is shown in Fig. 16. It is manifest that, for initial values 
of the chemical potential no % t above the singularity, the final renormalized value of \x lies very close to the VHS. 
These results are important to assure that the enhancement of the instabilities due to the divergent density of states 
does not rely on fine-tuning the Fermi level to the VHS, as the chemical potential tends to pin itself in a natural way 
to the singularity. 

1 

0.75 
0.5 
0.25 
(X 
-0.25 
-0.5 
-0.75 

0.25 0.5 0.75 1 1.25 1.5 1.75 2 
A 

FIG. 16. Scaling of the chemical potential as a function of the high-energy cutoff. The results correspond to the Hubbard 
coupling U = At. 



The integrals of Eq. (55) can be used now to find the solutions of Eqs. (10) and (11) displaying the superconducting 
instability. The form of the flow in the coupling constant space is shown in Fig. 17. In the case of bare repulsive 
interactions, cither the BCS couplings scale to zero for Vj > Vu, or there is an unstable flow giving rise to the 
superconducting instability when Vj < Vu- The latter instance is realized in lattice models which have a nearest- 
neighbor attractive interaction V besides the on-site U repulsive interaction. When V < 0, the bare coupling Vj = 
U + 4V is obviously smaller than the bare coupling Vu = U — 4V. We are however more interested in the case of the 
pure Hubbard model, in which the bare couplings lie in the diagonal of the first quadrant in Fig. 17. 

The couplings read directly from the hamiltonian of the Hubbard model correspond to the boundary between the 
regions of stable and unstable flow. This means that the slightest perturbation may drive the system to either of 
the two sides, which stresses the role played by the irrelevant operators under these conditions. There are actually 
perturbations that fade away when the theory is scaled to low energies, but that may be important because they may 
destabilize the flow in the BCS channel. 




15 




FIG. 17. Flow of the renormalized BCS couplings in the (Vi,Vu) plane. 

In the particular case of the Hubbard model, such irrelevant perturbations are given by the iteration of particle-hole 
diagrams of the type shown in Fig. 18. Apart from the particle-particle diagrams, these are the only corrections that 
arise from the bare couplings of the model, and they are not enhanced at low energies since the particle-hole bubbles 
do not have the appropriate kinematics to be of order ~ dA in the wilsonian approach [2]. 



ki -k t k! -k t 




Pt -Pi Pt -Pi 

(a) (b) 
FIG. 18. Particle-hole corrections to the BCS vertices in the Hubbard model. 

The iteration of the bubbles in Fig. 18 gives rise to antiscreening diagrams, i. e. to corrections that add to the 
bare repulsive interaction. We recall that the particle-hole bubble with total momentum about Q is enhanced with 
the factor d given after Eq. (8), while that with momentum about the origin is proportional to the factor c given 
after Eq. (3). As long as in the present paper we remain in the range t' < 0.276 t, we have that d is greater than c, 
and we face the instance in which the irrelevant perturbations make Vjj slightly larger than Vi at the beginning of 
the RG flow. 

We have solved the RG equations (10) and (11) taking as initial values for Vj and Vu the result of adding the 
ladder series built from the diagrams in Fig. 18, with a bare Hubbard coupling U = At. Moreover, in the resolution 
we have introduced the dependence of fi on A that arises from Eq. (55). This is one of the main accomplishments 
of our RG procedure, since the knowledge of how the VHS is approached is essential to regularize the effect of the 
divergent density of states. 

The results can be synthesized in the determination of the line at which the transition to the superconducting 
state takes place in the model. That is characterized by the energy at which the BCS couplings grow large or, more 
conveniently, by the point at which these couplings have a singularity. This depends on the initial position fio of the 
chemical potential, and it has been represented as a function of this variable in Fig. 19. 

We find that the BCS couplings diverge only for values of no in the range of attraction to the VHS, that is when 
the renormalized chemical potential is pinned to the singularity. There is an optimal value of no for which the scale of 
the transition reaches a maximum, as the chemical potential stays closer to the VHS during a greater part of the RG 
flow. For lower values of no, the scale of the instability decreases, as a consequence of the fact that the renormalized 
chemical potential is not precisely pinned then to the VHS. 
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FIG. 19. Plot of the energy scale of the superconducting instability (thin line) and of the transition to the spin instability 
phase (shaded region). 

We have also represented in Fig. 19 the energy at which the spin instability opens up, according to the estimate of 
Section V. We realize that this scale is always above the energy at which the singularity develops in the BCS channel, 
whenever the spin instability exists in the system. This happens for values of /iq higher than the optimal one. For 
lower values, the renormalized chemical potential deviates from the VHS by an amount even larger than the gap that 
would be due to the spin instability, so that this does not find the conditions to develop. We have then a picture 
in which the pairing instability exists alone for /iq below the optimal doping, but it is actually precluded above that 
level since the spin instability sets in before with the formation of a gap in the quasiparticle spectrum. 

We comment finally on the symmetry of the condensate wavefunction. The fact that the Umklapp interaction Vjj 
becomes increasingly repulsive when approaching the instability implies that the wavefunction must have opposite 
signs in the saddle-points A and B. As long as in the unstable flow we approach the asymptotic regime Vj = —Vu, 
the response function for the d-wave operator 

*+ T (k)*+ (-k) - *+ T (k)*+ ; (-k) + h.c. (59) 

develops a singularity at the frequency where the coupling Vj — Vu blows up. By the same token, it is easily seen that 
the response function for the s-wave operator does not display any divergence at low energies. Without the need of 
knowing precisely the shape of the gap, we may assure then that the symmetry of the order parameter is of d-wave 
type, with nodal lines at the bisectors of the four quadrants. This is in agreement with the results of more general 
analyses, which show that the symmetry of the order parameter can be ascertain from the topology of the Fermi line 
alone [12]. 

VII. CONCLUSIONS 

In this paper we have presented a study of the different phases of the system of electrons interacting near a Van 
Hove singularity, when the correlations at momentum Q = (jr, tt) prevail over those at zero momentum. In the context 
of a model with nearest-neighbor and next-to-nearest-neighbor hopping, this happens for < t' < 0.276 t, according 
to the comparison of the prefactors c and d that appear in Eqs. (3) and (8), respectively. We have applied a wilsonian 
RG approach following the same lines developed by Shankar in Ref [2] for the analysis of Fermi liquid theory. We 
have paid attention to the spin degrees of freedom when considering the different interactions, what has allowed us 
to discern the universality classes of the system. 

We have seen that, regarding the spin correlations, there is a universality class characterized by a spin instability 
in the low-energy theory, in opposition to the regime of couplings with smooth behavior of the correlators for the spin 
operators. In the case of the extended Hubbard model with on-site interaction U and nearest-neighbor interaction V, 
the spin instability arises for U > 0, irrespective of the value of V, and it is absent for U < 0. 

Several authors have previously considered the competition between the spin and the superconducting instabilitcs 
in the universality class corresponding to the divergent flow in the upper half-plane of Fig. 12 [15,17,25,26,29-31]. Our 
analysis has shed light into a number of features of the spin instability. We have seen that this takes place through the 
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condensation of particle- hole pairs with momentum Q. The fact that a macroscopic number of these pairs has been 
formed is what forces the quasiparticles to live out of the range A excited by the condensate in the neighborhood of 
the saddle-points. The Fermi line is destroyed in a region whose size is bounded by yJ\A\/t', in units of the inverse 
lattice spacing. 

This effect provides a paradigm for the disappearance of the Fermi surface of an electron system which differs 
from the understanding of such a phenomenon in Mott-Hubbard insulators. Those systems are supposed to be in 
a strong-coupling regime, in which the double occupancy of each lattice site is highly suppressed. In our case, we 
need otherwise to constrain the Fermi level near the VHS, departing sensibly from half-filling as if is increased. Most 
remarkably, the instability takes place no matter how small the bare couplings may be in the above picture. This is 
what ultimately allows to discern the symmetry breaking in the ground state within our RG approach. 

We have seen that two different behaviors arise also in the space of couplings for the BCS channel, starting from 
bare repulsive interactions Vi and Vjj- The (i-wave superconducting instability develops in models corresponding to 
the region with unstable flow in the upper half-plane of Fig. 17. This is the case of the extended Hubbard model with 
U > and attractive interaction V. When U > and the nearest-neighbor interaction is repulsive, the couplings 
in the BCS channel scale down to zero. The Hubbard model with just on-site interaction is placed at first sight on 
the boundary between the regions with stable and unstable behavior. We have shown that the model has irrelevant 
perturbations that drive the system towards the side with divergent RG flow. Since the departure from the limit 
behavior is weak, the superconducting instability is overshadowed by the spin instability, up to a point of optimal 
doping beyond which the latter is absent. 

The use of the wilsonian RG approach provides some advantages over other RG methods, the most important being 
the possibility of studying the renormalization of the chemical potential. Given the divergent behavior of the density 
of states at the VHS, it is clear that all the positions of the Fermi level cannot be equally stable. The scaling of 
the chemical potential can be obtained by letting it free to evolve and computing the shift from the integration of 
high-energy modes near the cutoff at each RG step. Following this procedure, we have seen that there is a range of 
attraction near the VHS where the chemical potential is renormalized down to the singularity. This guarantees the 
naturalness of the different instabilities since, rather than relying on the fine-tuning of the Fermi level, they arise from 
its precise pinning to the VHS in the low-energy effective theory. 



APPENDIX 



In this section we compute the imaginary part of some of the susceptibilities of the model. It turns out that the 
particle-hole susceptibility % p ? l (Q,w) and the particle-particle susceptibility Xpp(0,w) have a nontrivial imaginary 
part, while this vanishes for Xph(0,oj) and Xpp(Q> w ) at any finite frequency. 

In our model, the susceptibility Xph(Q> w ) is given by 

(Q u) = i [ <kjq ( - - (60) 

Xph ,UJ % J 2tt J (2tt) 2 lj + uj q - £A(q) +«e sgn(w + uj q ) uj q - e B (q) + «e sgn(w 9 ) 

where the energy cutoff is implicit in the integration over the momenta. According to the standard prescription, the 
imaginary part of (60) is given by 

Im Xph(Q^) = ~ 2n2 j Tj^T J (|^2 Sgn ( w + w «) s g n K) <K W + u i ~ £ a{<\)) 5{u q - £E?(q)) (61) 
In the limit of small t', Eq. (61) leads to a quantity which does not depend on the frequency. Taking w > 0, we have 

Im X ph(Q, u) = f d 2 q 5(l0 + 2t(q 2 x - q 2 y )) (62) 
1 [ q ° 1 

= ^J- qo dqx 7W+wW) (63) 

where qo = \J ^ - After a little of algebra, we obtain 



Im Xp/l (Q,c) = -log^- + ^-lj (64) 
c'/(87rf) (65) 
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We see therefore that the imaginary part is equal to tt/2 times the prefactor of log(A) in the real part of the 
susceptibility. It can be checked that the same relation holds between the real and the imaginary part of the particle- 
particle susceptibility x PP (0,w). 

Turning now to the susceptibility Xph(0,u>), we have 

f diOq f d 2 q 1 1 

Im Xvh{0, to) — Re / — — - I — — — - (66) 

J 2tt J (2-kY w + uj q - e A (q) + it sgn(o> + uj q ) uj q - e A {q) + ie sgn(Wg) 

= -2tt 2 j J (^yi sgn ^ + u i"> s s n K) S ( LJ + uj q - Mq)) 5 H - £ ^(q)) ( 67 ) 

= -±-5(u) J d\ (68) 

We see then that the imaginary part of the susceptibility is zero for any finite value of the frequency 

A result similar to (68) is obtained for the imaginary part of the susceptibility Xpp(Q) w )- I n this channel, the pole 
that arises after summing up the ladder series corresponds to the appearance of excited states in the spectrum. We 
conclude therefore that the breakdown of symmetry through a mechanism of condensation can only take place in the 
particle-particle channel at zero momentum and in the particle-hole channel at momentum Q, as stated in the text. 
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